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Abstract 

When a bubble of air rises to the top of a highly viscous liquid, it forms a 
dome-shaped protuberance on the free surface. Unlike a soap bubble, it bursts 
so slowly as to collapse under its own weight simultaneously, and folds into 
a striking wavy structure. This rippling effect occurs in fact for both elas- 
tic and viscous sheets, and a theory for its onset is formulated. The growth 
of the corrugation is governed by the competition between gravitational and 
bending (shearing) forces and is exhibited for a range of densities, stiffnesses 
(viscosities), and sizes ~ a result which arises less from dynamics than from 
geometry, suggesting a wide validity. A quantitative expression for the num- 
ber of ripples is presented, together with experimental results which are in 
agreement with the theoretical predictions. 



Every day, nature surprises us with structures and patterns of such beauty as to fill the 
scientist with wonder and the artist with envy. In the present paper, we address an instability 
which turns a hemispherical, smooth, liquid bubble into a striking wrinkled structure, first 
observed by Debregeas, de Gennes, and Brochard-Wyart In their experiment, 0.1 to 
10 cm^ of air injected into a highly viscous liquid (77 ~ 10'^ Pa-s) rises to the free surface, 
imprisoned in a hemispherical bubble of thickness t ~ 1 — 10 fim. If the latter is punctured 
at its apex by a needle, surface tension drives the rapid expansion of a circular opening. 
The retraction velocity soon (after about 10 — 30 ms) saturates to a constant, owing to 
the high viscous resistance. In the meantime, the air fiow through the hole equilibrates the 
pressure difference, allowing the bubble to collapse under its own weight. As it deflates, an 
instability appears: the fluid sheet folds into a wavy structure, with radial ripples that break 
the original axisymmetry. In the absence of a detailed theory, Ref. [|[] proposes a scaling 
estimate for the number of ripples n* ~ {figR^ / kY^'^ , where /i is the mass of the fllm per 
unit area, g the gravitational acceleration, R the radius of the hole, and K an effective 
bending rigidity of the sheet which was assumed to be elastic during the early stages of the 
rippling. 

The rippling results from the competition between compression, bending, and gravity. 
Each fluid element tends to fall under its own weight, but experiences a viscous resistance 
from its neighborhood. If the bubble were to collapse in a uniform, symmetric way, it 
would occupy a progressively reduced area, leading to an in-plane compression which would 
require forces that far exceed the scale set by gravity. Instead, the fllm deforms in a nearly 
inextensional fashion by undergoing pure bending. Equivalently, for a given (gravitational) 
force, the relative time scale associated with stretching is much larger than that for bending, 
and the surface therefore corrugates over short times, before eventually relaxing into a 
uniform, thicker membrane. 

This instability is reminiscent of buckling phenomena [Q], originally studied in the con- 
text of elastic rods but occurring also in the creeping flows of viscous liquid fllaments (a 
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striking everyday example being the coiling of a stream of honey when it reaches a toast 
For an elastic rod, buckling occurs at the longest possible wavelength in order to 
minimize the bending energy. In the bubble problem, however, gravity plays a distinctive 
role in determining the configuration. For a given amplitude, bending still favors large scale 
deformations, while gravitational energy is minimized for an almost flat sheet with as many 
tiny ripples as possible; the optimal wavelength results from a compromise between the 
two. Such an argument, however, does not fully characterize the effect. Unlike the above 
examples, here the system under consideration is a curved two-dimensional sheet, and the 
associated geometry constrains the rippling both qualitatively and quantitatively, as we shall 
see. 

The instability occurs both an elastic (solid) film and a viscous (liquid) one. The elastic 
case corresponds to a convex curved shell with an apical hole of radius R that is allowed to 
collapse under its own weight. In the viscous case, an additional complication arises because 
the radius of the hole changes during the rippling. Following a short initial transient, the hole 
grows steadily at a rate v ^ a/r] resulting from the balance of surface tension {a ~ 20 mN/m) 
and viscous stress. It thus takes a time r ~ f]t/'^ the opening radius to increase by t. 
During this time, the hquid acquires a velocity V gr due to gravity, larger than v by 
a factor V/v ^ 10^. Even if the liquid is viscoelastic, so that the retraction velocity is 
enhanced by a factor R/t (~ 10 — 10^) the hole radius remains essentially constant 

while the instability occurs (c/. Fig. 1). We may therefore treat the hole radius i? as a given 
parameter in the theory P|. 

Although the bubble has the geometry of a sphere before collapsing, it is quite flattened 
by the time the ripples appear (Fig. 1(B)). For simplicity, we consider the unperturbed 
configuration to be a shallow cone, of slope a <^ 1, described by its height above the 
surface, 

h = a{ro-r), (1) 

where r is the cylindrical radial coordinate and tq the radius of the base. Any deformation 
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of h introduced by the rippling may be written, without loss of generality, as 

h + Sh = a{ro-r) -6a (r) + ^ [6(3^^^ (r) cos (nO) + 6(3^^^ (r) sin (nO) } , (2) 

where 6 is the azimuthal angle. The perturbation 6a represents a uniform (n independent) 
flattening accompanying the growth of ripples of amplitude 6(3^^ , and a crucial step consists 
in understanding their form and interdependence. In the case of a thin elastic (viscous) sheet, 
the two primary modes of deformation are in-plane stretching (shearing) and out-of-plane 
bending. A generic deformation of an elastic cone (made of a material with Young modulus 
Y), of amplitude C on a scale i, requires stretching forces (per unit surface) of order Yt^^/i"^, 
but significantly smaller stretching forces (per unit surface) of order Yt^C,'^/l'^ 0, so that 
for a given external force, here gravity, inextensional deformations are greatly preferred [Q. 
In the case of a highly viscous sheet, forces arise from velocity gradients, thus introducing a 
dynamical element in the problem. However, their dependence on t and £ (essentially due to 
the variation of the strain across the film) is similar, so that inextensional deformations are 
again largely favored ii t tq/ti* . (This condition is satisfied if the selected number of 

ripples n* is small compared to 10'^, which is the case as we shall see below.) Equivalently, for 
a given loading, the time scale corresponding to bending is smaller than that for stretching 
by a factor (t//)^ |]9|,p^. Thus, at the onset of the instability, perturbations of the cone must 



preserve its metric. This requirement translates into the constraints ITT] 



<5/3« (r) = 6pf X r + (3) 
(where 6(3!^^ and 6(3^^' are constants) and 

Aa6a (r) = ^ | - l) ^f^f" - ^o) + n'5f3f' (---))• (4) 

In the following, we elucidate the elastic (solid) case before extending our treatment to 
the viscous (liquid) case. The energy functional of a perturbed elastic cone reads 



E \h + ih] = / i (surface) X (gravitational potential energy + bending potential energy) 

J cone 

= ^2 rf^Vl + Nih + 6h)f {h + 6h) + ^ 'I . (5) 



where K = Yt^ / 12 (1 — z/^) is the rigidity and v the Poisson ratio. Only bending elastic 
energy appears in E, since we have confined ourselves to the class of inextensible deforma- 
tions. 

If the elastic cone is attached to the plane on which it rests so that 5h (r = tq) = 0, 
Eq. (|) yields 5^'^'^' = -6/3^^ro. On substituting Eqs. into Eq. (g) we then obtain, to 

lowest order in the perturbation, 

5E = E[h + 6h]- E [h] 



7T 

2' 



= E ^En- (6) 



n>l 

7"^/'^ = [aK j iigf^"^ is an intrinsic length scale arising from the competition between gravity 
and bending elasticity |]T2 . 



Each mode contributes an amount bE^ to the change in energy, and rippling occurs 
if bEn < for some integer. In general, 6En < for a range of different n's; the most 
negative variation corresponds to the maximally growing perturbation, and thus sets the 
wavelength of the instability. The formulation also yields a "threshold condition" ip [r^/ R) < 
'-fR'^ip {tq/ R) for the occurrence of rippling. This condition involves the three independent 
quantities 7, tq, and R, and may be translated into three corresponding statements, (i) 
Rippling is suppressed if 7 < 7c (ro, -R) = R'^ip/ip, i.e. if the cone is too light or too 
rigid, (a) Similarly, no rippling occurs if the hole, or equivalently the cone, is too small, 
'"o < ''"oc {l,R)- Azimuthal continuity requires the wavelength of the deformation to be at 
most of order tq, resulting in a forbidding bending cost if tq becomes small compared to the 
intrinsic (energetically determined) scale 7"^/^. {iii) Finally, the threshold depends, quite 
unexpectedly, on the ratio tq/R. The dependence of the symmetric (n = 0) mode on the 
radial coordinate r is different from that of the rippling {n 7^ 0) modes, so that the high 
elastic cost cannot be justified by gravitational gain anymore if the hole is reduced beyond 
a critical size. Minimizing 6E in Eq. (10) yields the selected number of ripples as 
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where Int{x} is the integer closest to x. This relation improves on the estimate of Ref. [|l|, 
where the authors consider the short time elastic behavior, and establishes its domain of 
validity. 

For an elastic (solid) sheet, the rippling phenomenon is of an essentially static nature; 
upon increasing, say, the mass of the sheet, the equilibrium configuration is shifted from 
symmetric to rippled. Approaching the problem from a dynamical perspective by considering 
the elastic forces and torques rather than the corresponding energies results in an evolution 
equation 7r/iP(ro,/2) • £ /dt^ = -5E/6f3^^ for each mode. Here P is a polynomial 

function independent of n, so that the energetically optimal mode, with number n*, is indeed 
the fastest growing one. In the case of a viscous liquid, the effect is intrinsically dynamical: 
bending occurs only on short times, while the equilibrium configuration is ultimately reached 
by a slow thickening. Nevertheless, the motion of a viscous film satisfies a formulation close 
to that of an elastic sheet, as can be shown by integrating the Stokes equation through 
the thickness [|10|,|13[. Indeed, it is easy to see that bending results from a torque ^i^^2-^ x 
(i(curvature) /(it analogous to an elastic torque K x (curvature), so that a highly viscous 
film may be described by an effective bending modulus Ki = rjt^/ST (z/ = 1/2 for an 
incompressible medium) , where r is a time scale associated with the falling velocity. Thus, all 
the conclusions of the stability analysis for the elastic cone, and in particular the expression 
for the number of ripples (Eq. (|^)), may be transposed to the case of the bubble modulo a 
certain time scale related to the gravity-induced velocity of the fluid. Comparing the nascent 
ripples' amplitude to the film thickness yields an estimate of this time scale as (t/gY^'^ [14 



In order to check our results against experiment, we visualized the bursting of silicone 
oil bubbles. Once the bubble is punctured with a sharp needle, its evolution is followed 
using a high-speed camera capable of recording up to 1000 frames per second. The resulting 
video is then analyzed to determine the radius tq of the bubble, the hole size R at which 
the ripples are first observed, along with the number of ripples n*. Since the hole expands 
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very fast at first, R is much larger than Rc by the time the bubble begins to collapse. To 
compare the experiments with the theory, in which R enters as a parameter, the latter is 
measured at the onset of the instability for each given size of the bubble. The quantitative 
measurements are compared to the theoretical predictions for the dependence of n* on the 
bubble size on Fig. 2. On a more qualitative level, the experiments show a suppression of 
the instability for small bubbles, in agreement with the threshold conditions above. 

We conclude with a discussion of possible refinements of the theory and their relation to 
the geometric nature of the problem. A more complete theory would incorporate a (flattened) 
hemisphere as the initial condition, rather than a cone. Also, due to the progressive drainage 
of the liquid, the thickness t acquires a dependence on r (and time). This in turn implies 
non-uniform rigidity K (r) and mass fi (r), leading to functions /, g, (p, ip, and P of a more 
complicated form. On a more fundamental level, all these aspects should be addressed in 
terms of the coupled hydrodynamics of the slow viscous (liquid) flow and the rapid air flow; 
the corresponding calculations will be reported elsewhere ||13|. Yet, the strong geometrical 
constraints involved in the problem are suggestive of the robustness of the results. The 
question we have answered is akin to that of applying a curved surface unto a fiat one in the 
most economical way, a problem which has taxed cartographers for many centuries and lies 
at the birth of differential geometry. It is also somewhat of an inverse counterpart to the 
problem of fitting a flat sheet to a three dimensional landscape, which has been studied in 



various contexts [p!5|-[r7|, and is an issue that still vexes fashion designers. The relevance of 
the geometrical constraints is manifest, for example, in the strong dependence of the rippling 
on the size of the opening, which is closely related to a well-known theorem due to Gauss 
181 , Jellett [|19[, and others, according to which (loosely put) a closed surface cannot be bent 



without being stretched, while an open surface can be bent inextensionally. Similarly, we 
find that a smaller hole implies a relatively stiffer bubble, and hampers the rippling. While 
the detailed form of the functions if and arise from the physical constraints and dynamics 
imposed by the forces and various boundary conditions, the essence is in the geometry. 
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FIG. 1. Stroboscopic images of a collapsing liquid bubble of size is tq = 1cm and thickness 
t ^ 100/im. The silicone oil has viscosity t] = lO^Pa ■ s, surface tension cr = 21mN/m, and 
mass density 0.98g/cm'^. (A) 30ms after the film is punctured by a sharp needle, the bubble 
shows a retracting hole of radius R = 1.4mm at its apex, but no ripples yet. (B) 30ms 
later, the bubble loses its axisymmetric shape. The radius of the hole remains essentially 
constant, at i? = 1.6mm, while the ripples grow. The inset displays a schematic side view 
of the essentially conical deflating bubble at the onset of the instability, with the important 
quantities involved in the phenomenon. The extreme shallowness allows for a perturbative 
treatment in the slope a of the cone. 

FIG. 2. Plot of the number of ripples n* as a function of the bubble radius tq, comparing 
the experimental measures (points) with the theoretical predictions (solid lines). These data 
were gathered using silicone oil of viscosity t] = 600Pa-s and bubbles of thickness t ^ 30/im. 
The errors in the measurement of tq arise from meniscus effects which are more important 
in smaller bubbles. The bursting time elapsed up to rippling is measured to be of order 
of 1 to 5 times (t/g)^^"^, consistent with our proposed mechanism for the formation of the 
corrugation. For each experimental realization, the ratio tq/R was measured at the onset of 
the instability, and the corresponding dependence of R on vq was used to get a theoretical 
curve n* = n* (ro). The thinner line displays the prediction for an elastic sheet attached to 
the plane on which it rests. The liquid, however, can in no way be clamped, and one has to 
relax the boundary conditions at the base. This leads to a vanishing of the unprimed modes 
(c/. Eq. (^), which are unfavorable in terms of both graviational and viscous forces; the 
fastest growing primed modes lead to the behavior represented by the thicker line. The latter 
is plotted here for a slope a ^ 3°(~ O.OSrad) of the cone, consistent with our perturbative 
treatment, and in agreement with direct observation. The dashed line represents the best 
fit of the scaling form n* ~ (figR'^/K)^^'^ where R is chosen as the relevant length 
scale. If R is replaced by ro, the above expression for n* may be closely fitted (up to an 
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overall multiplicative factor) to our predicted curve, showing that the size of the bubble is 
the dominant length scale within the present experimental range and conditions. This is 
consistent with the relaxed boundary conditions, which allow the ripples to be significant 
close to the outer edge of the bubble (c/. also Fig. l.(B)). In this way, the ripples trade 
a bulk gain in gravitational and bending stresses against a cost in stretching in a thin rim 
close to the outer edge. The increased thickness of the liquid film close to the base further 
emphasizes this effect, as it reduces the difference in magnitude between a typical stretching 
and a typical bending stress. The inset shows a top view of the fully developed ripples, from 
which n* is measured. 
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This figure "bubble_after.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0008175vl 



This figure "bubble_before.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0008175vl 
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This figure "bubble_top.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0008175vl 



